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ABSTRACT 

We examine the long term (^10 years) 7-ray variability of blazars observed by EGRET and 
Fermi and find that for six sources the average flux varied by more than an order of magnitude. 
For two of these sources (PKS 0208-512 and PKS 0528+134), there were extensive observa- 
tions (at various observing periods) by EGRET. Hence these dramatic variations are not due 
to a single short time-scale flare, but reflect long term changes in the average flux. Over the 
last twenty years, these two sources were also the target of several X-ray observatories (e.g. 
ROSAT, ASCA, RXTE, BeppoSAX, Chandra, Suzaku, XMM -Newton and Swift). While the 
ratios of the average 7-ray fluxes between EGRET and Fermi observations are 22.9 ± 1.9 
and 12.6 ± 1.5, their estimated 2 — 10 keV X-ray flux do not show such dramatic variations. 
The X-ray emission from such flat spectrum radio quasars (FSRQs) are believed to be due 
to synchrotron self Compton, while 7-rays originate from inverse Comptonization of external 
soft photons from an accretion disk and/or broad line region. We argue that in this scenario, 
the only explanation for the uncorrelated variability is that there was an order of magnitude 
decrease in the external soft photons, while the jet parameters remained more or less constant. 
This result indicates that perhaps the accretion and jet processes are not tightly coupled in 
these sources. 

Key words: galaxies: active — galaxies: jets — quasars: individual: PKS 0208-512 — 
quasars: individual: PKS 0528+134 — gamma-rays: galaxies — X-rays: galaxies 



1 INTRODUCTION 

Blazars are a subclass of active galactic nuclei (AGN) whose 
bolometric luminosity is dominated by 7-rays which sometimes 
extends to TeV energie s. According to the unification scenario 
dUrry &Padovanilll995l) . they are jet dominated sources having a 
small jet to line-of-sight angle. They show high variability through- 
out the electromagnetic spectrum, particularly in 7-rays which are 
believed to be emitted close to the base of the jet. However, the 
nature and e mission mechanisms of the jet are still not w ell un- 
derstood (e.g. lFerrarifl998t lB6ttcherJl2007l; lMarscherll2009l) . 7-ray 
observations of these sources can provide a better understanding 
of fundamental issues relating to the energetics, origin of the jet 
emission and possible connection to the underlying accretion disk. 

COS-B sate llite detected the first extragalactic 7-ray 
source (3C 273; iBignami et al] Il98ll) . The first all sky sur- 
vey in 7-ray was c o nducted by E G RET onboard CGRO 
dKanbach et all 1 19881 : iHartman et ail 1 19991) . EGRET de- 
tected 271 discrete sources above 100 MeV and ~100 
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of the identified sou r ces of EGRET catalogue are blazars 
( Hartman etalj 1999; Sowards-E mmerd. Romani. & Michelsonl 
120031 : ISowards-Emmerd et alj 120041) . Based on their optical line 
emission properties, blazars can be divided into two sub classes: 
flat spectrum radio quasars (FSR Qs) and BL Lacs with th e FSRQs 
being more luminous (e.g. lKrolilll999l : lKembhavilll999l) . EGRET 
detected significantly higher number of FSRQs than BL Lacs. The 
Fermi 7-ray space telescop e (Fermi), l aunched in 2008, has ~ 30 
times better sensitivity dAtwood et all l2009h and sky coverage 
than EGRET, and hence it has detected ~ 1900 sources from its 
more than two y ears of observations (Second Fermi catalogue: 
iNolan et alj|2012l) . Most of the identified Fermi sources are AGNs 
and the clean sample of 2nd LA T-AGN catalogue contains 886 
sources dAckermann et al.l 1201 lb . Majority of these AGNs are 
blazars - 395 BL Lacs, 310 FSRQs, and 157 candidate blazars of 
unknown type. 

There have been several a ttempts to mod el spectral energy 
distributions (SEDs) of blazars dBoettch er 2004 and the references 
therein). For example , model ling the SED of the well known FSRQ 
3C 279 ISikora et al.l d200lh . concluded that X-rays and 7-rays 
might be produced co-spatially by electrons of similar energies. 
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They considered both synchrotron self Compton (SSC) and ex- 
ternal Compton (EC) emission processes in their model and their 
results indicated that while X-rays are mainly produced via the 
SSC, 7-rays are dominated by EC process. This has been con- 
firmed by using more detailed models and for other FSRQs (e.g., 



Table 1. Details of highly variable blazars 



Abdo et alj boiOal; iPalma et al J I20TTI ; ISahavanathan & Godambel 
2012l ; lMukheriee etal.ll 1999b . 



Blazars show variability in 7-rays over different time scales 
ranging fro m flares as short as ~ 1 hour to flares extending for days 
or months (jHartman et al.lll999t iMattox et alii 19971; iHartman etal] 



l200ll:lBoettcherll2004l) . Blazars are also variable on still longer 
timescales of years or decades. The most striking evidence for 
this is that a fraction of the EGRET detected blazars are still not 
detected by the more sensitive Fermi telescope. Moreover, some 
EGRET blazars (e.g., PKS 0208-512) were found to be mostly in 
the low state during Fermi observations. While there have been 
several studies regarding the rapid variability and flaring activity 
of bl azars and their possible connection to varying jet param eters 
(e.g. iBottcher et alj|2003l ; ISahavanathan & Godambel 120121) . it is 
not clear whether the very long term (~ 10 years) variations are 
similar or if they are intrinsically different phenomena. 

In this work we study the long term variability (~ 10 year) of 
FSRQ that were monitored by EGRET during its operation (1991- 
1997) by comparing their average EGRET flux with those mea- 
sured recently by Fermi i.e. after 2007. From 1997 to 2007, there 
was no 7-ray observatory available to study these sources at GeV 
energies (7-ray satellite AGILE was launched in 2007, but no long 
term monitoring has been reported). There were pointed X-ray ob- 
servations of these sources by different satellites (e.g., ROSAT, 
RXTE-PCA, BeppoSAX, SWIFT, etc.). Our motivation is to ex- 
amine any correlation between the 7-rays and X-rays, which in 
turn would provide important clues to the main driver of the long 
term variability. We concentrate on the 7-ray and X-ray variability, 
since SED modelling indicates that although they both arise from 
the same region (indeed the same electron distribution), their radia- 
tive mechanisms are different and hence their response to variation 
in th e intrinsic jet proper ties will be different. 

lAbdo et all l l2010bh and IPalma etal] ( 1201 ll) studied the rela- 
tively short term variability o f blazars uti l ising qu asi-simultaneous 
multiwaveband observations. I Abdo et a l. (2010b) carried out quasi 
simultaneous multiwaveband observations of 5 7-ray loud blazars 
between 2008 October and 2009 January. Using the result from 
SED modelling they argue that the difference between the low and 
high state of luminous blazars is due to the different kinetic power 
of the jet which is most likely related to the varying bulk Lorentz 
factor of outflow within the blazar emission zone. Palma et al. 
j201ll) reported that PKS 0528+134 did not show any significant 
variability in 7-ray band on 1 — 2 weeks timescales. They found 
moderate variability in X-rays ( | ^f- | ~ 50%) and most radio fre- 
quencies (]^-| < 20%) on timescales of 1 — 2 weeks. In optical 
bands the variability is found to be of upto AR < 1 on timescale 
of several hours. However, the physics driving the long term vari- 
ability may be different than the short term one. 



2 OBSERVATIONS 

From the 3rd EGRET point source catalogue dHartman et al.ll999h . 
blazars have been identified using multi-waveleng t h info rmation 



jHartman etal] 1 19991 ; IMattox. Hartman. & Reimerl 1200 ll). Using 
a diff erent identification technique Sowards-Emmerd et alj d2003l 
2004) have obtained a more complete census of plausible blazar 



Source Name 


h EGRET 


T, a 

r Fermi 


R a 


PKS 0208-512 


86.7 ±4.6 


3.79 ± 0.24 


22.9 ± 1.9 


3EGJ 1614+3424 


10.1 ± 1.5 


0.49 ±0.11 


20.6 ±5.5 


PKS 1622-29 


40.1 ±3.1 


2.40 ±0.26 


16.7 ±2.2 


PKS 0336-01 


21.4 ±5.0 


1.35 ±0.17 


15.9 ±4.2 


PKS 0528+134 


32.4 ± 1.2 


2.57 ±0.29 


12.6 ± 1.5 


PKS 1454-354 


1.1 ± 0.3 


11.9 ±0.42 


0.09 ± 0.02 



a 7-ray (> 1 GeV) flux during EGRET and Fermi observations in 10 9 
photon cm~ 2 s~ 1 .R is the ratio of the EGRET to Fermi flux. 
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Figure 1. EGRET lightcurve (> 100 MeV) of PKS 0208-512 (top panel) 
and PKS 0528+134 (bottom panel) averaged over each viewing period. The 
dash line indicates the average P1234 flux value. MJD 48200 and MJD 
50000 correspond to 1990 November 5 and 1995 October 10 respectively. 



counterparts. We use their source list augmented with source iden- 
tifications from 3rd EGRET catalogue in regions of the sky that 
were not covered bv lSowards-Emmerd et al. : (2003. 2004). We fur- 
ther consider only FSRQs that have a cumulative source detection 
significance > 4cr for sources above the Galactic plane (\9\ > 10°) 
and > 5cr for those in the Galactic plane (—10° ^ 6 ^ 10° ). The 
7-ray fluxes averaged over all the cycles of EGRET observations 
(PI 234) were used for our calculations. The PI 234 flux has been 
calculated from the counts map over four cycles. Our final source 
list contains 56 FSRQs out of which 38 are also present in the 2nd 
Fermi catalogue. To be compatible with the Fermi data, we con- 
sider only EGRET 7-ray fluxes above 1 GeV, which are calculated 
from the average spectral index and > 100 MeV flux quoted in the 
3rd EGRET catalogue. 

We find 5 FSRQs which are brighter by an order of magnitude 
or more during EGRET era compared to the average flux observed 
from Fermi data while there is one source that is fainter by a factor 
of ~ 10 (Table 1). 

We focus first on two of these highly variable six sources (PKS 
0208-512 and PKS 0528+134) that have been detected in several 
EGRET viewing cycles. These two sources were also observed in 
X-rays by different satellites during EGRET to Fermi era. Figure 
Q] top and bottom panels show EGRET lightcurves of PKS 0208- 
512 and PKS 0528+134 respectively. The dash line indicates the 
P1234 average flux value. The figure shows that there was extensive 
coverage of these sources by EGRET and the average values are not 
dominated by a single, short time-scale flare. 
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Table 2. Summary of X-ray observations of PKS 0208-512 and PKS 0528+134 



Satellite Date Energy band (keV) Flux a Flux2_iofceV 6 Ref. c 



PKS 0208-512 



ROSAT 


1990-1991 


0.1 - 


2.4 


1.72 ±0.36 


0.98 ±0.21 


1 


ASCA 


1995 


0.5 - 


- 10 


9.49 


6.76 


2 


ASCA 


1996 


2 - 


10 


5.77 


5.77 


3 


BeppoSax 


2001 


2 - 


10 


4.7 


4.7 


4.5 


Chandra 


2001 


2 - 


10 


2.1 d 


2.1 d 


6 


SWIFT 


2005 


0.5 - 


- 10 


3.4 


1.9 


7,8 


SWIFT 


2008 


0.5 - 


- 10 


2.7 


1.6 


7,8 


Suzaku 


2008 


2 - 


10 


1.37 ±0.06 


1.37 ±0.06 


7 



PKS 0528+134 



ROSAT 


1990-1991 


0.1 - 


- 2.4 


3.13 ±0.63 


3.99 ±0.80 


1,9,10 


ASCA 


1994 


2 - 


10 


o c+0.8 
°*°-0.4 


q c+0.8 


11 


ASCA 


1995 


2 - 


10 


8.98 


8.98 


3 


RXTE-PCA 


1996 


2 - 


10 


10.62 ± 0.70 e 


10.62 ±0.70 e 




BeppoSax 


1997 


2 - 


10 


2.5 


2.5 


5,12 


RXTE-PCA 


1999 


2 - 


10 


7.00±0.58 e 


7.00±0.58 e 




Suzaku 


2008 


2 - 


10 


2.8 


2.8 


13 


RXTE-PCA 


2009 


2 - 


10 


< 7.85 e 


< 7.85 e 




XMM-Newton 


2009 


0.2 


-2 


n 9 o+0.03 


1.35 ±0.11 


13 


SWIFT 


2009 


0.2 - 


- 10 


4.0 ±0.6 


3.0 ±0.5 


14 



"Reported X-ray flux in 10 12 erg cm 



in the energy band given in column 3. 'Using the fluxes and spectral parameters observed in different energy 



bands given in referenc es the X-ray fluxes a r e conv ert ed in 2- 10 keV en erg y band in 10 12 erg cm 2 s 1 unit. c Reference s: 1: 
iBrinkmann. Siebert. & BollerNl994l) 2:|Reeves & Turner! fcOOOl) 3:IUeda et"al]j200lh 4:lTayecchio et alji l 2002l) 5 : Donato. Sambruna. & Glioz3 fe005l) , 6: 
"lAbdo et alifeOlObl) 8:IZhang et all fepiOl) 9:IZhang et alji'l994l) 10:lMukheriee et aljfl996h ll:ISambruna et all jl997|) 12: 



Schwartz et al 
lGrrisellrnTet*alT 



fe006l> 7: 
1999) 13: 



Palma et all 1201 ll) 14: lGhisellini, Tavecchio. & Gh irlanda (2009) Flux has been calculated from the rest-frame 2 — 10 keV 



luminosity given in the reference. e Averaged over the observations shown in Fig|2] 
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Figure 2. RXTE-PCA lightcurve of PKS 0528+134. The left-hand branch 
of points of the top panel represents the observations during 1996 August 
and September, and the middle and right-hand branch represent observa- 
tions during 1999 January and May, respectively. The bottom panel repre- 
sents observations during 2009. The horizontal dash line indicates the aver- 
age flux over the subset of data. 



PKS 0528+134 is one of the most luminous blazars de- 
tected by EGRE T and showed se veral flaring episodes from 1990 
to 1997 iMukh eriee et alj 1 1997b . It was also detected by both 
OSSE and COMPTEL onboard CGRO jMcNaron-Brown et alj 
ll995tlCollmaretalJ l997) and has been ob served simultaneously at 
other wavelengths ( Mukherie e"e7aT]ll999l and references therein). 
PKS 0208-512 shows the maximum flux variations with an aver- 



age EGRET flux that is 22.9 ±1.9 times higher than the average 
Fermi flux. Over the last twenty years, both sources were observed 
by several X-ray satellites and Table [2] gives their X-ray fluxes as 
reported. These different instruments measured X-ray fluxes in dif- 
ferent energy bands (column 3 of Table|2}- Using the reported spec- 
tral parameters, we converted the observed flux in different energy 
bands into flux in the 2-10 keV band. Specifically, we use the re- 
ported column density and spectral index to model the spectrum in 
XSPEC and extracted the 2-10 keV flux (column 5 of Tabled}. 

Pointed observation by RXTE Proportional counter array 
(PCA) of PKS 0528+134 were undertaken in 1996, 1999 and 2009. 
We used the spectra available in the standard product of the PCA 
observations, modelled with a power-law and Galactic absorption 
to obtain a 2 — 10 keV integrated flux. Figure|2]shows the flux vari- 
ations across different observations. During 2009, only upper limits 
were found probably because only one of the Proportional Counter 
Unit was working during this time. During these observations the 
flux did not vary significantly and hence we report an average value 
for three periods of observations in Table |2] 

Figures [3] compares the 7-ray and X-ray fluxes at different 
times for the two sources. EGRET 7-ray fluxes are averaged over 
each observing period (~ 2 weeks) and the Fermi flux is aver- 
aged over the first two years of observation. X-ray fluxes are from 
pointed observations. While there is considerable scatter in the X- 
ray fluxes, there is no evidence for a dramatic decrease in X-ray 
flux, especially for PKS 0528+134. For PKS 0208-512, there is a 
hint of a decrease in the X-ray flux if one neglects the early ROSAT 
point. However, this decrease is at most a factor of 3, which is sig- 
nificantly less than the variation in the 7-ray flux of a factor of 23. 
Thus, we conclude that the large variability observed in 7-rays is 
not seen in X-rays. 
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Table 3. Summary of X-ray observations of other four FSRQs 



Satellite 


Date 


Energy band (keV) 


Flux" 


Flux2_iofc e v b 


Ref. c 


PKS 1622-29 












ROSAT 


1990-1991 


0.1 - 2.4 


0.32 ±0.08 


0.31 ±0.08 


1,2 


RYTF-PfA 


1997 


2 — 10 




ii ar. _|_ n ok 




RXTE-PCA 


2006 


2-10 


14.49 ±0.25 


14.49 ±0.25 


- 


PKS 0336-01 












Einstein 


1979-1980 


0.15 - 3.5 


0.48 ±0.20 


0.68 ±0.28 


3 


ROSAT 


1990-1991 


0.1 - 2.4 


0.98 ±0.29 


0.73 ±0.22 


1 


PKS 1454-354 












ROSAT 


1990-91 


0.1 - 2.4 


0.510 


0.290 


1.4 


SWIFT 


2008 January 


2-10 


0.6 ±0.12 


0.6 ±0.12 


5 


SWIFT 


2008 September 


2-10 


1.8 ±0.4 


1.8 ±0.4 


5 


3EGJ 1614+3424 












ROSAT 


1990-91 


0.1 - 2.4 


0.401 


0.228 


1,4 



"Reported X-ray flux in 10 12 erg cm 2 s 1 in the energy band given in column 3. 6 Using fluxes and spectral parameters observed i n different energy 
bands given in r eferenc es X-ray fluxes are deri ve d for a common 2-10 keV energy b an d in 10~ 12 erg cm~ 2 s — 1 unit. c References: l:lBrinkmann et all 
il994l) 2: lMattoxetai]h997h 3: l0wen. Helfand. & Spangled J 198 lh 4: lMassaro etllU2009l) 5: lGhiselhni et al.H2009l) 



For the other four sources with high 7-ray flux variability, only 
sparse X-ray observations are available. Nevertheless we report all 
available X-ray fluxes in Table[3] For PKS 1622-29, we have anal- 
ysed existing RXTE PCA observations of 1997 and 2000. As com- 
pared to the early ROSAT observation, X-ray flux has increased 
in contrast to 7-ray emission which decreased. However, the small 
number of observations do not permit a reliable conclusion. For the 
rest, there are no reported X-ray observations during both EGRET 
and Fermi era. 

Statistical errors on the measured flux have been reported for 
some of the observations. However, due to the uncertainties in the 
relative normalisations between different instruments and in the ex- 
trapolation to the 2-10 keV band, systematic errors on flux could be 
higher, but not expected to be larger than a factor of 2. 

For two of the blazars, PKS 0528+134 and PKS 1622-29, 
the RXTE All Sky Monitor (ASM) lightcurves are available in 
the RXTE website (http://xte.mit.edu/asmlc/One-Day.htmll. How- 
ever, the long term time averaged count rates for these sources PKS 
0528+134 and PKS 1622-29 are not significa ntly higher than th e 
systematic positive bias expected in ASM data dChitnis et al.l2009h . 
Thus ASM did not detect these sources. 

It is to be noted that X-ray coverage of these sources are quite 
sparse and hence, the actual short term variability can be under- 
sampled. However, our main aim is to search for any long term vari- 
ability which is less affected by the under-sampling of light curves. 
Also, for PKS 0528+134, quasi-simultaneous gamma-ray and X- 
ray observations have been carried out during almost all EGRET 
observations. It was observed that apart from a single huge flare, 
the variation in both X-ray and gamma-ray emission during these 
times are less than a factor of 3. 



3 DISCUSSION 

In leptonic jet models, non-thermal electrons in the jet produce 
synchrotron emission. These electrons upscatter the synchrotron 
photons (SSC) and also photons coming out from accretion disk 
and/or broad line region (BLR) (EC) via the inverse Comp- 



ton process. These emission processes contribute to the SED of 
blazars and modelling suggests that 7-ray emissions of FSRQs 
are dominated by EC proc ess while the X-rays are primarily 
due to SSC (e.g., 3C 279 : iPian et alj 1 19991 : lAbdo et al.1 l2010al: 
Saha vanathan & Godamba I2012L PKS 0528+134: IPalmaetalJ 
201 lh . On the other hand, SEDs of BL Lac sources are better fit- 
ted by synchrotron self Compton (S SC) emission (e.g., Mrk 501: 
IPetrv et al.ll200ol lAbdo et allkoiOah . Recent pro gress in theoreti - 
cal models for blazar emission is summarised bv lBoettcherlteOlfjl) . 

Specifically, Mukherj ee et al.1 dl999h carried out a multiwave- 
band study of PKS 0528+134 during EGRET observations and 
SED modelling showed that during the 7-ray high state, the 7- 
ray emission was d ue to externa l Com pton (EC) emission. Using 
Fermi observations, IPalma et"al, I J201 1) also carried out a coordi- 
nated multi-wavelength campaign and similarly concluded that as 
expected for FSRQs, the 7-ray emission is mainly due to EC pro- 
cess while SSC process is the primary X-ray emission mechanism. 

In this scenario, we investigate variability in which jet or 
environment parameters can give rise to a strong variation in 
the EC component (i.e. 7-rays) while not affecting the SSC one 
(i.e. X-rays). Variation of the magnetic field density is ruled out 
since that would affect SSC and not the EC emission. Similarly, 
a change in the number density of electron n would affect SSC 
more (oc n 2 ) than the EC emission (oc n). Since, blazar emis- 
sion is highly beamed, observed flux is a sensitive function of 
the Doppler beaming factor (5). For a discrete jet, synchrotron 
self Compton emission is proportional to 8^ 3+a \ whereas, exter- 
nal Compton flux is proportional to fi( 4 + 2a \ where a is the en- 
ergy spectral index. The average value of 7-ray photon spectral 
index (r = a + 1) derived from Fermi and EGR ET data sets 
are 2.47 ± 0.19 and 2.34 ± 0.15 respectively dAbdo et alJkoiOct 
Bhattacharva, Sreeku mar. & Mukherjee 2009) and hence variation 
in a cannot be the cause of the 7-ray intensity variation. For the 
average value of a = 1.4, the Doppler factor S would need to de- 
crease by ~ 1.4 to cause a factor of 10 decrease in the EC flux. 
However, this would have led to a decrease in the SSC component 
as well by a factor of ~ 5, which was not observed. Thus, varia- 
tion in any of the intrinsic jet parameters such as number density, 
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day (MJD) 

Figure 3. X-ray and 7-ray fluxes of PKS 0208-512 (top panel) and PKS 
0528+134 (bottom panel) as observed by different satellites. The stars rep- 
resent X-ray observations from different instruments from EGRET to Fermi 
era. RXTE fluxes are averaged over observations shown in Fig |2] It is ev- 
ident that average X-ray flux did not decrease by more than a factor of 2 
from EGRET to Fermi era. The triangle represents the 7-ray flux (> 1 
GeV) averaged over first two year of Fermi observations. The boxes repre- 
sent EGRET fluxes (> 1 GeV) averaged over each observing period (~ 2 
weeks). The scale for gamma-ray fluxes are given in the right side of the 
figure. MJD 48000 and MJD 56000 correspond to 1990 April 19 and 2012 
March 14 respectively. 



magnetic field, spectral index or Doppler factor cannot explain the 
uncorrelated variation in X-ray/7-ray emission. 

The EC emission is proportional to the external photon flux 
from accretion disk or BLR region, while SSC emission is inde- 
pendent. If the external photon flux was a factor of ~ 10 more dur- 
ing EGRET observation time, then it would naturally explain why 
7-rays changed by a factor of ten while the X-ray flux remained 
nearly constant. This has important implications on the relationship 
between the jet and the accretion process. If the flux from the disk 
or BLR, which is connected to the accretion process can change by 
an order of magnitude but the average intrinsic properties of the jet 
are not affected, it implies very weak or no coupling between the 
accretion and jet emission processes. 

For PKS 0528+134, the X-ray observations showed no signif- 
icant variation i.e. the variation was less than a factor of 2. If one 
considers the extreme conservative limit, wherein the X-ray emis- 
sion did decrease by a factor of 2 during the Fermi era, it implies 
that the Doppler factor did decrease by a factor ~ 1.17. However, 



this decrease in Doppler factor alone cannot explain the observed 
large 7-ray flux variation. An additional decrease in the external 
photon flux by a factor of ~ 5 is required. Thus in this extreme 
case, a factor of 5 change in the accretion process caused only a 
modest ~ 1.17 factor change in the Doppler factor, implying that 
the accretion and jet processes are not tightly coupled. 

If the external photons arising from the accretion disk or the 
BLR region did indeed change significantly over time, one can ob- 
serve a significant decrease in the optical broad line flux during the 
Fermi era compared to EGRET. Such a decrease, would be a di- 
rect evidence for the model proposed here. However, in order to 
obtain an averaged line flux long term spectroscopic monitoring of 
these sources in optical is required. While, ther e were several pho- 
tometric observations of PKS 0528+134 ( e.g. IVillata et alii 19971; 
lRaiterietal.il 19981; iMukheriee etai1ll999h . optical s pectroscopic 
obser vations have been carried out only in recent years dPalma et alj 
l201ll) . 

Apart from synchrotron process, the photometric observation 
will have contributions from accretion disk, BLR region and also 
from the host galaxies. So the optical photometric flux is not al- 
ways domi nated by synchrotron and/or synchrot r on self Compton 
proce sses {whiting, Maiew ski. & Webstedl2003l IMukheriee et alj 
Il999l) . PKS 0528+134 was observed in optical during few EGRET 
observing periods and the variation is within a factor of ~ 3. This 
source was also observed in optical during the n ight of 2001 Au- 
gust 29 jWhiting et al.l2003h . |Palma et aljfeoi lh also observed the 
source in 2009 and the variation between these observations (BVR 
band) is a factor of ~ 2. 

There were only few optical observations of PKS 0208-512. 
This source was obse rved in optical wav ebands during the nights 
of 2001 August 29 dWhiting et ai1l2003l) and 2008 December 14 
dAbdo et al.l2010bh . So, it is not possible to have a meaningful long 
term average of optical emission from this object. 

PKS 1622-29 was monitored in optical (mainly in R band) 
for more than a decade by the program for Extragalactic Astron- 
omy (PEGA) dOsterman Meyer et al. 2008). Though the source 
showed strong variation in s horter time scale (Figure 1 and 2 of 
lOsterman Meyer e t al. 2008]), the long term average variation is 
within a factor of ~ 3. However, the number of pointed observa- 
tions is too few to derive a reliable long term average flux. More 
importantly, a huge 7-ray flare of this source was observed with 
EGRET and hence, the calculated EGRET average flux may be 
overestimated. Therefore, we have mainly considered PKS 0208- 
512 and PKS 0528+134 for our analysis. Unfortunately, as per our 
knowledge, there are no reported optical spectra of these blazars 
during both EGRET and Fermi times. 

Radio emission probably originates in a different emission re- 
gion than that of X-ray and 7-ray. Hence, the high frequency ra- 
dio flux might be delayed by a few months to that of X-rays and 
7-rays. Since the main objective of this work is to study the long 
term behaviour of blazars (~ 10 year), due to the longer cooling 
time scales at radio wavelengths, the high-frequency radio power 
could be a good representative for the long time-averaged jet power. 
PKS 0528+134 has been observed at a wavelength of 850^im from 
1997 April to 2004 Decem ber using SCUBA camera at the Ja mes 
Clerk Maxwell Telescope jjenness, Robson. & Stevensll2010h . A 
total of 129 observations were carried out over 62 nights. The 
so urce did not show large variations during this period ( - 



3). lLonsdale. Doeleman. & Phillips! 11998) conducted a survey of 
bright extragalactic radio sources with 3 mm (100 GHz) VLBI. 
They reported that the 3 mm fluxes are 2.9 Jy and 2.5 Jy as mea- 
sured by Kitt Peak and Haystack telescopes and MPIR 100 me- 
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ter telescope near Bonn and Onsala telescopes respectively. This 
source has been studied in the observed frequency range 80 — 267 
GHz from 2007 February to 2008 December with IRAM Plateau 
de Bure Interferometer (PdBI). Reported average fluxes in 1.3 mm 
band, 2 mm band and 3 mm ban d are (2 ± 1.1) Jy , (1.7 ± 0.5) 
Jy and (3.6 ± 1.3) Jy respectively jMaurv et alfe oiO). Clearly, the 
flux variation in 3 mm band between th ese two observations is less 
than a factor of 2. lAgudo et alj ( bold) observed PKS 0528+134 
during 2005 July at 86 GHz and the reported flux is (2.0 ± 0.10) 
Jy. This source was also monitored at 86 GHz during 20 01 Oc- 
tober and the reported flux is 2.02 Jy fcxe et alj 12008). This 



ent times and no large flux variation is noticed (e.g., IPalma et al 


2011 




Massardi et al. 2009^ Gu, Cao. & Jians 


2009; Lanyietal 


2010 


; 


WirenetalJ 1992; lAbdo et al. 2010a; 


Petrov et al.l 120071; 


Lister & Homan 2005; Kellermann et al. 2004; 


Cara & Listej2008 



etc.). 

The measured 150 GHz flux of PKS 0208-512 is 1268 ± 
86 mjy utilising the observation with Arcminute Cosmology 
Bolo meter Array Rece iver (ACBAR) during 2002 observing sea- 
sons dKuo et alj[2007h . WMAP has detected PKS 0208-512 and 
no significant variation is noticed among one y ear averaged, 
three year averaged and five year averaged fl uxes dBennett et alj 
l2003l;lHinshaw et alj|2007l; I Wright et alj2009l) . Using the five year 
WMAP data, IChen & Wrightl d2009h reported that observed 41 
GHz, 61 GHz and 94 GHz fluxes of this source are (2.7 ± 0.2) 
Jy, (2.7 ± 0.2) Jy and (1.8 ± 0.4) Jy respectively. PKS 0208- 
512 has been observed at other high radio frequ encies at different 
times and no large flux variation is noticed (e.g., Ricci et al]|2006l ; 



iMassardi et alj2008l ; lMurphv et al.l2010l ; lTingav et al.ll2003f ; etc.). 

The results obtained in this work need to be confirmed for 
other highly variable 7-ray sources through reliable long term X- 
ray monitoring data. It would also be interesting to study long term 
behaviour of such sources at other wavelengths which may provide 
additional clues on the strength of coupling between the jet and 
accretion process. 
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